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1. Introduction 

It is well established that at high baryon density the combination of asymptotic freedom and 
the existence of attractive channels in the color interaction between the quarks lying in the large 
Fermi surface come together to promote the formation of quark-quark pairs, which in turn break the 
color gauge symmetry giving rise to the phenomenon of color superconductivity. At densities much 
higher than the masses of the u, d, and s quarks, one can assume the three quarks as massless and 
the favored state results to be the so-called Color-Flavor-Locking (CFL) phase characterized 
by a spin zero diquark condensate antisymmetric in both color and flavor. 

The conditions of extremely large density and very low temperature required for color super- 
conductivity cannot be recreated in Earth's labs. Fortunately, nature provides us with a laboratory to 
probe color superconductivity, the cores of celestial compact objects. These compact stars typically 
have very large magnetic fields. Neutron stars can have magnetic fields as large as B ~ 10 12 — 10 14 
G in their surfaces, while in magnetars they are in the range B ~ 10 14 — 10 15 G, and perhaps as high 
as 10 16 G [§] (for a recent review of magnetic fields in dense stars see Even though we do not 
know yet of any suitable mechanism to produce more intense fields, the virial theorem [Q] allows 
the field magnitude to reach values as large as 10 18 — 10 19 G. If quark stars are self-bound rather 
than gravitational-bound objects, the upper limit that has been obtained by comparing the magnetic 
and gravitational energies, could go even higher. 

A natural question to ask is: What is the effect, if any, of the huge star's magnetic field in the 
color superconducting core? A complete answer to this question would require a rather involved 
study of quark matter at the intermediate range of densities proper of neutron stars, where the 
strange quark mass cannot be ignored, with the additional complication of an extra parameter, the 
magnetic field. However, as a first, more tractable approach to this question, one can ignore the 
strange quark mass effects and look for the consequences of an external magnetic field on the 
superconducting phase, assuming that the quark matter is formed by three massless flavors. This 
was the strategy followed in our recent paper [jfj|], whose main results will be described in what 
follows. 

In this talk I will show the way a magnetic field affects the pairing structure and hence its 
symmetry, ultimately producing a different superconducting phase that we have called Magnetic 
Color-Flavor-Locking (MCFL) phase. 

In a conventional superconductor, since Cooper pairs are electrically charged, the electromag- 
netic gauge invariance is spontaneously broken, thus the photon acquires a Meissner mass that can 
screen a weak magnetic field, the phenomenon of Meissner effect. In spin-zero color superconduc- 
tivity, although the color condensate has non-zero electric charge, there is a linear combination of 
the photon and a gluon that remains massless This new field plays the role of the "in-medium" 
photon in the color superconductor, so the propagation of light in the color superconductor is dif- 
ferent from that in an electric superconductor. 

Because of the long-range "rotated" electromagnetic field, a spin-zero color superconductor 
may be penetrated by a rotated magnetic field B. Although a few works |Q] had previously ad- 
dressed the problem of the interaction of an external magnetic field with dense quark matter, none 
of these studies considered the modification produced by the field on the gap itself. However, as 
we have recently shown the gap structure gets modified due to the penetrating field. To un- 



2 



Dense Quark Matter in a Magnetic Field 



Vivian de la Incera 



derstand this, notice that, although the condensate is <2-neutral, some of the quarks participating 
in the pairing are Q-charged and hence can couple to the background field, which in turn affects 
the gap equations through the Green functions of these ^-charged quarks. Due to the coupling 
of the charged quarks with the external field, the color-flavor space is augmented by the Q-charge 
color-flavor operator, and consequently the CFL order parameter splits in new independent pieces 
giving rise to a new phase, the MCFL phase. 

2. The MCFL Gap Structure 

The linear combination of the photon and a gluon G^ that behaves as a long-range field in 
the spin-zero color superconductor is given by [[I], R], 

A M = cos 0A M -sin0 G®, (2.1) 

while the orthogonal combination = sindA^ + cos$G^ is massive. In the CFL phase the mix- 
ing angle 6 is sufficiently small (sin0 ~ e/g ~ 1/40). Thus, the penetrating field in the color 
superconductor is mostly formed by the photon with only a small gluon admixture. 

The unbroken £7(1) group corresponding to the long-range rotated photon (i.e. the t/(l) e m .) 
is generated, in flavor-color space, by 2 = 2x1 — lx Q, where Q is the electromagnetic charge 
generator. We use the conventions Q = — As/v^, where As is the 8th Gell-Mann matrix. Thus 
our flavor-space ordering is (s,d,u). In the 9-dimensional flavor-color representation that we will 
use in this paper (the color indexes we are using are (l,2,3)=(b,g,r)), the Q charges of the different 
quarks, in units of e = ecos 6, are 
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(2.2) 



In the presence of an external rotated magnetic field the kinetic part of the quarks' Lagrangian 
density must be rewritten using the covariant derivative 

l%arks = V{^t)V , (2-3) 

with 

n M = id M + efiA M . (2.4) 

where 

Q = Q.+ - Q. . (2.5) 
is the rotated charge operator. The charge projectors 

Q.+ = diag (0,0,0,0,0,0,1,1,0), (2.6) 

=diag(0, 0,1, 0,0, 1,0,0,0) , (2.7) 
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and 



obey the algebra 



Q = diag(l, 1,0, 1,1,0,0,0,1) 



(2.8) 



0.^0.^ = 8^0.^, 17,77' = 0,+, 



(2.9) 



O + &+ + ^ = 1 • 



(2.10) 



The rotated magnetic field naturally separates the quark fields according to their Q charge. 
The fermion field in the 9x9 representation used above, \j/ T = (si,S2,s^,di,d2,d^,ui,U2,UT,), can 
then be written as the sum of three fields with zero, positive and negative rotated electromagnetic 
charges, 

V = V(o) + V(+) + V(-) , (2- 11 ) 
where the (0)-, (+/— )-charged fields can be respectively written in terms of the charge projectors 
as 

y {0) = OqW , = O+y , Y(-)=^-¥- (2-12) 

A strong magnetic field affects the flavor symmetries of QCD, as different quark flavors have 
different electromagnetic charges. For three light quark flavors, only the subgroup of SU (3)l x 
SU (3)r that commutes with Q, the electromagnetic generator, is a symmetry of the theory. Simi- 
larly, in the CFL phase a strong B field should affect the symmetries of the theory, as Q does not 
commute with the whole locked SU (3) group. Based on the above considerations, and imposing 
that the condensate should retain the highest degree of symmetry, we proposed [§] the following 
ansatz for the gap structure in the presence of a magnetic field 



/ 2a; 











Af 







A s - A A 




















A A +A 5 
As-A A 



Aa 
















2A 5 




Af + Af 



Af 















•A2 















Af- 















•A2 



A« + Af\ 




A^ + Af 



2A C 



(2.13) 

We call the reader's attention to the fact that despite the Q-neutrality of all the condensates, 
they can be composed either by neutral or by charged quarks. Condensates formed by Q-charged 
quarks feel the field directly through the minimal coupling of the background field B with the 
quarks in the pair. A subset of the condensates formed by <2-neutral quarks, can feel the presence 
of the field via tree-level vertices that couple them to charged quarks. The gaps A^ , s are anti- 
symmetric/symmetric combinations of condensates composed by charged quarks and condensates 
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formed by this kind of neutral quarks. The gaps A A , as well as As, A s and A s , on the other hand, are 
antisymmetric and symmetric combinations of condensates formed by neutral quarks that do not 
belong to the above subset. The only way the field can affect them is through the system of highly 
non-linear coupled gap equations. At zero field the CFL gap matrix is recovered since in that case 
A^ = A A and Af = A s = A' s = A" s . 

Although the symmetry of the problem allows for four independent symmetric gaps, the con- 
densates A s and A^ are only due to subleading color symmetric interactions, and as explained in 
the previous paragraph, they are formed by neutral quarks that are not coupled to charged quarks, 
so they belong to the same class as A5. Therefore, there is no reason to expect that they will differ 
much from A5. Hence, in a first approach to the problem, we will consider A5 ~ A^ ~ A^. 

The order parameter ( |2.13 ) implies the following symmetry breaking pattern: SU (3) co i OT x 
SU(2) L x SU(2) R x U(1) B x U^(\) A x C/(l) e . m . -> SU (2) color+L+s x U(l) em , The U^(l) A 
symmetry is connected with the current which is an anomaly-free linear combination of s,d and u 
axial currents [|0. The locked SU (2) corresponds to the maximal unbroken symmetry, and as such 
it maximizes the condensation energy. Notice that it commutes with the rotated electromagnetic 
group t7(l) em . 

The counting of broken generators, after taking into account the Anderson-Higgs mechanism, 
tells us that there are only five Nambu-Goldstone bosons. One is associated to the breaking of the 
baryon symmetry; three Goldstone bosons are associated to the breaking of SU (2) A , and another 
one associated to the breaking of U^\\) A . All the Nambu-Goldstone bosons are <2-neutral. The 
number and properties of the lightest particles in the MCFL have implications for the low-energy 
physics of the phase. Since in her talk Cristina Manuel will address the low-energy physics of the 
MCFL phase, I will not extend on this topic in mine. 

3. Effective Action in the Presence of a Magnetic Field 

Let us construct the effective action of the system in the presence of a magnetic field. With this 
aim, we will use a Nambu-Jona-Lasinio (NJL) four-fermion interaction abstracted from one-gluon 
exchange Although this simplified treatment disregards the effect of the fi-field on the gluon 
dynamics and assumes the same NJL couplings for both the situation with and without magnetic 
field, it keeps the main attributes of the theory, thereby providing the correct qualitative physics. 

We start from the mean-field effective action 



h{w,w)= {jW{o)(x)[gj q)q \ (x,y)V(o)(y) + Y(+){x)[G( +)0 ] (*,y)V(+)(y) 



+i/ H (x)[G+ )0 ] l (x,y)Y(-){y) + v mc (x)[G m) } 1 (x,y)V(o)cCy) 
+¥(+)cW[ G (+)o] _ H^^)^(+)c(3')+¥(-)cW[ G (l)or 1 ^,3 ; )^(-)c(3')] 

+\ [?(o)cW A+ ^y) V(o) 60 + h.c] + - [y {+)c (x)& + (x,y) ¥(-) (y) 

+W(-)c( x ) A+ (^y)V(+)(y)+h-c.}} , (3.i) 

where the external magnetic field has been explicitly introduced through minimal coupling with the 
Q— charged fermions. The presence of the field is also taken into account in the diquark condensate 



A + = 75A, whose color-flavor structure is given by Eq.(2.13). 



5 



Dense Quark Matter in a Magnetic Field 



Vivian de la Incera 



In (3.1) symbols in parentheses indicate neutral (0), positive (+) or negative (— ) Q— charged 



quarks. Supra-indexes + or — in the propagators indicate, as it is customary, whether it is the 
inverse propagator of a field or conjugated field respectively. Then, for example, [G^q] -1 corre- 
sponds to the bare inverse propagator of positively charged conjugate fields, and so on. The explicit 
expressions of the inverse propagators are 

[G% )Q r\x,y) = [ifd, l -m±llf]8\x-y) , (3.2) 
[Gf +)0 ]-\x,y) = [ifU^-m±^]8\x-y) , (3.3) 



[Gf_ )0 )-\x,y) = [ifn^-m±iif}S\x-y) , (3.4) 



with 



Ylf ] =id ll ±eA ll . (3.5) 

Transforming the field-dependent quark propagators to momentum space can be performed 
with the use of the Ritus' method, originally developed for charged fermions [Q] and recently 



extended to charged vector fields [ 1C]. In Ritus' approach the diagonalization in momentum space 
of charged fermion Green's functions in the presence of a background magnetic field is carried out 
using the eigenf unction matrices E p {x). These are the wave functions of the asymptotic states of 
charged fermions in a uniform magnetic field and play the role in the magnetized medium of the 
usual plane-wave (Fourier) functions e' px at zero field. Below we introduce the basic properties of 
this transformation. 

The transformation functions Eq (x) for positively (+), and negatively (— ) charged fermion 
fields are obtained as the solutions of the field dependent eigenvalue equation 



(Jl^-Y)E i t\x)=E ( t\x){y-p^) , (3.6) 



with given by 



and 



with eigenfunctions 



p^ = (p Q ,0,±^2\eB\k,p 3 ), (3.7) 
Ejf\x)=^E£ ) (x)8(o), (3.8) 



E$\x) = ^e-^+^+P^D^p^) , (3.9) 
where D„ ±) (P(±)) w& the parabolic cylinder functions with argument P( ± j defined by 



P{±] = j2\eB\(x l ±p 2 /eB), (3.10) 
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and index nr±\ given by 



2\eB\ z 



H±) 



0,1,2, 



(3.11) 



k = 0,1,2,3,... is the Landau level, and a is the spin projection that can take values ±1 only. 
Notice that in the lowest Landau level, k = 0, only particles with one of the two spin projections, 
namely, a = 1 for positively charged particles, are allowed. The normalization constant ^Ki {±) is 



jY n{±) = {A%\eB\Y/Jn {±) 



In (3.8) the spin matrices 8(o) are defined as 



8(a) =diag(5 ff i,5 ff _i,5 ff i,5 ff _i), 
and satisfy the following relations 



o = ±\ , 



5(±) t = 5(±) , 8(±)8(±) = 8(±) , 5(±)5(=f) = 0, 



f8(±) = 8(±)f, r L 8(±) = 8(T)7 L 



In Eq. ( p.!5[ ) the notation y" = (y°, y 3 ) and y- 1 = (7 1 , y 2 ) was used. 
The functions are complete 

£ f dp dp 2 dp 3 E^\x)Ef\y) = (2n) 4 8^(x-y) , 

and orthonormal, 

/ Ef\x)E<t\x) = (27i) 4 A k 8 kk/ 8(po-p'o)8( P2 -p' 2 )8(p 3 -p' 3 ) 
with the (4 x 4) matrix A k given by 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



A* 



8(0 = sgn[eB}) 
I 



for 
for 



k = 0, 
k>0. 



(3.18) 



The matrix structure A,t was recently introduced in Ref. [11]. It had been previously omitted 
in the orthonormal condition of the E p (x) functions given in Refs. [Kj 12]. Nevertheless, 
it should be underlined that this matrix only appears in the zero Landau level contribution, and 
consequently it enters as an irrelevant multiplicative factor in the Schwinger-Dyson equations in 
the lowest Landau level approximation. Thus, all the results obtained in the works [|[ |l0[ 12] 



remain valid. In Eqs. ( |3.16| )-( 3T7| ) we introduced the notation E p (x) 



UE^ix^yo. 



Under the E p (x) functions, positively (V(+))> negatively charged fields transform ac- 

cording to 

V(±)M =E / dp dp 2 dp i E p ± \x)Y(±)(p) , (3.19) 
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W(±)(x) = £ / dp dp 2 dp 3 y {±) (p)E p ±) (x) . 
k J 

One can show that 



and 



[Y^(-)^±^ Q )-m]E ( p \x)=E p ; ±ju5 M0 )-wi] 
The conjugate fields transform according to, 



,(-). 



V(+)cW =L / dpodp 2 dp 3 E ( p \x)y (+)c (p), 
k J 

V(-)cto =L / dp dp2dp 3 E p + \x)Y(-)c(.P) ■ 



(3.20) 

(3.21) 
(3.22) 

(3.23) 
(3.24) 



After transforming to momentum space one can introduce Nambu-Gorkov fermion fields of 
different Q charges. They are the <2-neutral Gorkov field 



the positive 



and the negative one 



(o) 



*<+) = 



*(-) = 



V(0) 
V(0)C , 

V(+) 

V(-)c, 

V(-) 

V(+)c 



(3.25) 



(3.26) 



(3.27) 



Using them, the Nambu-Gorkov effective action in the presence of a constant magnetic field B 
can be written as 



+ 



1 / (p) ^ + 1 / (W^-) {p) ' (3 " 28) 



2 7 (2 



(2^) 4 
1 f d 4 p 



where 



'(0)0 
A 



(0) 



(3.29) 



(0) 



r (0)0J 



Hp) 



[G+ 



(+)oJ 



~\P) 



"(+) 

{ G M »r\p) 



(3.30) 
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(\GU\-\P) AT, 



'(_)0J WJ "(-) 
A(l) [G H0 ]-'(/» 



(3.31) 



with 



A+ +) =£1_A+£1 + , (3.32) 
Atj=n+A + n_, (3.33) 

A+ 0) = ^oA+ilo, (3.34) 

Notice that to form the positive (negative) Nambu-Gorkov field we used the positive (negative) 
fermion field and the charge conjugate of the negative (positive) field. This is done so that the ro- 
tated charge of the up and down components in a given Nambu-Gorkov field be the same. This way 
to form the Nambu-Gorkov fields is mandated by what kind of field enters in a given condensate 
term, which in turn is related to the neutrality of the fermion condensate (y c y) with respect to the 
rotated Q-charge. 

In momentum space the bare inverse propagator for the neutral field is 

[Gfo )o r l (p) = lMPn±^o)-m] , (3.35) 

where the momentum is the usual p = (po,Pi,P2,P3) of the case with no background field. 
For positively and negatively charged fields the bare inverse propagators are 

[Gf+XjrV) = M/i +) ±Mfyo) - m] , (3.36) 



and 



{ G Uo\^(P) = fetfi^o)-" 1 ] ( 3 - 37 ) 



respectively. 

4. Gap Solutions 

The main question we would like to address now is: Can we find a region of magnetic fields 
where the gaps A^ and A^, (or Aj and A|), differ enough from each other that the system is not in 
the CFL phase anymore, but in the MCFL phase? To explore the possible answer to this question 



we need to solve the gap equations derived from the Nambu-Gorkov effective action (3.28). 
In coordinate space the QCD gap equation reads 

A+( X ,y) = i^Xl f S21 (x,y)f hDf v {x,y) , (4.1) 

where S2i(jc,;y) is the off-diagonal part of the Nambu-Gorkov fermion propagator in coordinate 
space and, for simplicity, we have omitted explicit color and flavor indices in the gap and fermion 
propagator. Here Z>^ is the gluon propagator. 
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In a NJL model the gap equation can be obtained from Eq. (11) simply by substituting the 
gluon propagator by 



,8 AB 8^(x-y). 



(4.2) 



1 

The NJL model is characterized by a coupling constant g and an ultraviolet cutoff A. The 
ultraviolet cutoff should be much larger than any of the energy scales of the system, typically the 
chemical potential. In the presence of a magnetic field we should also assume that A is larger than 
the magnetic energy \fiB. In other studies of color superconductivity within the NJL model, the 
values of g and A are chosen to match some QCD vacuum properties, thus hoping to get in such 
a way correct approximated quantitative results of the gaps. We will follow the same philosophy 
here, noticing however that this completely ignores the effect of the magnetic field on the gluon 
dynamics. 



To solve the gap equation ( |4.1[ ) for the whole range of magnetic-field strengths we need to use 
numerical methods. We have found, however, a situation where an analytical solution is possible. 
This corresponds to the case eB > \i 2 /2. Taking into account that the leading contribution to the 
gap solution comes from quark energies near the Fermi level, it follows that for fields in this range 
only the LLL (Z = 0) contributes. 

Using the approximation A^ 3> A|,Aa, and Aa 3> As, the gap equations decouple and the 
equation for A^ is 

Aa ~ 3A 2 " ^ 



A2 



+ 



dq 



A 



g 2 eB 
3A 2 J- A (2k) 2 



A 



(4.3) 



where the first/second term in the r.h.s. of Eq.(fO|) corresponds to the contribution of Q-neutral/charged 
quark propagators, respectively. For the last one, we dropped all Landau levels but the lowest, as 
we are interested in the leading term. 



The solution of Eq. (4.3) reads 



A| ~ 2a/5A4 exp 



3 A 2 71 



2„2 



g z (H 2 + eB 



(4.4) 



with 8 = A — il. It can be compared with the antisymmetric CFL gap [13] 

3AV 



A CFL 



2\J 8\x exp 



2g 2 /i 2 

In this approximation the remaining gap equations read 
\B S 2 f d 3 q A r 



(4.5) 



6A2 ^ ( 2 *) 3 J{q-nf + 2{Aiy 



+ 



g 2 eB 



dq 



6A2 J - A (27r)2 J( q -n? + W 



(4.6) 
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g 2 f d 3 q (Yl A A 

Aa 



7 

9 I{q-Hf+2(^ A Y' 



+ AA ), (4-7) 



and 



<i 3 g / Aa 



18A 2 7 A (2,)3V^7— ^ 
-At) 2 + 2(Af)2 



(4.8) 



We express below the solution of these gap equations as ratios over the CFL antisymmetric 
and symmetric gaps 

Aa 1 / 36 21 1 3 \ 

exp - — + ———— + — , (4.9) 



a cfl 2 (7/34) 17* nx(l+y) 2* 

where x = g 2 \i 2 / 'A 2 n 2 , and y = eB/pL 2 , and 

Af _ A B A [3 9 v-1 



+ ^^—^ . (4-10) 



a cfl a cfl ^ 4 2xln2v+l 
A s A A 3 / 4 \ 
A^ ~ A™^ V~TVy) ' ( } 

Note that our analytic solutions are only valid at strong magnetic fields. The lower value 
eB ~ jU 2 / 2 corresponds to ~ (0.8 - 1.1) • 10 18 G, for jj, ~ 350-400 MeV. For fields of this 
order and larger the A^ gap is larger than A^ FL at the same density values. Nevertheless, we have 
estimated (see the details in Cristina Manuel's talk in the proceedings) that the separation between 
CFL and MCFL will take place already at fields ~ 10 16 G. 

All the gaps feel the presence of the external magnetic field. The effect of the magnetic field 
in A^ is to increase the density of states, which enters in the argument of the exponential as typical 



of a BCS solution. The density of states appearing in (4.4) is just the sum of those of neutral and 
charged particles participating in the given gap equation (for each Landau level, the density of 
states around the Fermi surface for a charged quark is eB/2n 2 ). 

All the <2-charged quarks have common gap A^. Hence, the densities of the charged quarks 
are all equal. As two of these quarks have positive Q charge, while the other two have it negative, 
the <2 neutrality of the medium is guaranteed without having to introduce any electron density. 

5. Conclusions 

In this paper, we have shown that a magnetic field leads to the formation of a new color-flavor 
locking phase, characterized by a smaller vector symmetry than the CFL phase. The essential 
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role of the penetrating magnetic field is to modify the density of states of charged quarks on the 
Fermi surface. To better understand the relevance of this new phase in astrophysics we need to 
explore the region of moderately strong magnetic fields eB < ji 2 /2, which requires to carry out 
a numerical study of the gap equations including the effect of higher Landau levels. Because the 
total density of states around the Fermi surface for charged particles does not vary monotonically 
with the number of Landau levels, we still expect to find a meaningful splitting of the gaps at these 
fields and therefore a qualitative separation between the CFL and MCFL phases. 
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